Engineering Entanglement: The Fast- Approach Phase Gate 
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Optimal-control techniques and a fast-approach scheme are used to implement a collisional control 
phase gate in a model of cold atoms in an optical lattice, significantly reducing the gate time as 
compared to adiabatic evolution while maintaining high fidelity. New objective functionals are given 
for which optimal paths are obtained for evolution that yields a control-phase gate up to single- 
atom Rabi shifts. Furthermore, the fast-approach procedure is used to design a path to significantly 
increase the fidelity of non-adiabatic transport in a recent experiment. Also, the entanglement power 
of phase gates is quantified. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.65.-w, 03.65. Ca 



I. INTRODUCTION 



Quantum information processing with cold atoms in 
optical lattices 0, 0, IE El IE IE an d m micro-traps 
DUE rely on the ability to entangle nearest neigh- 
bor atoms in an efficient controlled way. The effective 
trapping potential in optical lattices can have the form 
of a double well where the distance between minima, the 
well separation height, the frequency of each well, etc., 
can each be manipulated by controlling the lasers, using 
different combinations of laser frequencies and bias fields 
[To| . Ideally each atom is initially in one of the optical 
wells, and single qubits are registered into each atom's 
state. The atom's internal states (e.g.-hyperfine levels) 
|E0| or motional states in the trap [EEl can be used as 
a computational basis, \n), where n = 0, 1. A two-qubit 
gate is implemented by bringing the atoms together and 
letting their wave functions overlap [E EL IE IE IE IE IHI • 
During this overlap a phase that depends on the two- 
qubits state, \m, n) = \m)®\n), is accumulated because of 
the atom-atom (molecular) interaction. Two-qubit phase 
gates based on this scheme were suggested 0, |j, Q , ana- 
lyzed |E IE Sll, and demonstrated jllj . 

To achieve efficient computation and to be faster than 
decoherence processes, it is desirable to generate gates 
which operate as fast as possible 0, 0. Cirac and 
Zoller discuss the slow two-qubit collisional gate as 
one among two serious obstacles (the other is decoher- 
ence) to quantum computing with atoms in optical poten- 
tials. In designing a fast collisional gate one is faced with 
the problem of leakage outside the computational sub- 
space due to rapid switching of the control parameters. 
A theoretical question with immediate practical impact 
on the feasibility of quantum computation emerges: how 
fast can such a gate be? Previous theoretical work used 
adiabatic evolution to ensure the fidelity of the gate. In a 



recent experimental implementation [Tll j. the time scale 
for motion, 40 /xs, was chosen "to avoid any vibrational 
excitations" . 

Designing a fast two-qubit collisional gate is the pur- 
pose of this paper. Wc propose a fast-approach phase 
gate and use optimal-control to implement it. We first 
define the two-qubit <f> phase gate and explain its impor- 
tance. Next we review the adiabatic realization of such 
gates in timc-dependent optical potentials. Then we sug- 
gest the fast-approach scheme where adiabaticity is not 
required. Optimal control is applied and the nature of 
the resulting dynamical path is analyzed. Section [n] de- 
scribes two-qubit phase gates, Sec. 11111 sets out a simple 
model of a two-qubit phase gate, and Sec. lIVI presents the 
optimal control scheme for optimizing the gate thereby 
implementing a fast control phase gate. In Sec. [V] we 
present numerical results of the optimization, Sec. IVII 
applies the fast-approach technique to improve the gate 
experimentally demonstrated by and Sec. IVIll con- 
cludcs the paper. Appendix 1X1 shows that the phase <f> of 
a two-qubit phase gate uniquely the entanglement power 
of the gate, and Appendix IB1 describes how to implement 
a non-degenerate double well potential in an optical lat- 
tice. 



II. TWO-QUBIT GATE 

The two-qubit </> phase gate is designed to entangle two 
interacting atoms (qubits) to a desired degree: 

P(<j))\m, n) = exp {id mn ) \m,n) , 

#oo + #n - #oi - 0io = 4> (mod 2tt) . (1) 

This family of gates includes the control-phase gate: 



CP(0)|m, n) = exp (imncj)) \m, n) 



(2) 



'Professor Bilha Segev died on March 17 of complications arising 
in her battle with cancer. We lost a very dear colleague; a colleague 
who shared her enthusiasm for science and for life with us. She will 
be sorely missed. 



Any P(0) gate can be combined with unitary single- 
qubits Rabi shifts, exp(ia m ) exp(i/?„), to create the 
CP(<p) gate, where 6 mn + a m + /?« = mn(j) (mod27r), while 
#oo + #11 — #oi — #io is invariant under these shifts @ . The 
phase <p has an intrinsic physical feature in that it param- 
eterizes uniquely the entanglement power of the gate (see 
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Appendix A), and thus is intimately connected with the 
coupling strength of the two atoms during evolution. The 
control phase gate with phase (j> = it can further combine 
with single qubit operators to form the control-not gate, 
CN|m, n) = |m, n 771), where © denotes addition mod- 
ulo 2. Each of these two-qubit gates can be combined 
with a set of generators for single qubit gates to form a 
universal set for quantum computation [l5j. In practice, 
a physical system may evolve more naturally to a gate in 
P(-7r) other than CN or CP(7r). Therefore, in designing a 
gate, it is better to aim less restrictively for any one of 
the equivalent P(cf>) gates. 

The basic idea of realizing phase gates with a two par- 
ticle system in an external potential is: The external 
potential initially localizes the particles far enough apart 
so that they may be considered independent. The exter- 
nal potential then changes in time so that wave function 
overlap gives rise to correlations due to particle-particle 
interaction. The external potential is finally restored to 
its initial shape, so that the two particles no longer in- 
teract, but are now in a new correlated state. 



III. THE MODEL 

As in Refs. [H, IE IS EH 113 we focus on a two-qubit 
collisional gate where nearest-neighbor interaction is ma- 
nipulated by a time-dependent potential. A simple model 
for this process is a time-dependent Hamiltonian with a 
double well potential, whose minima are separated by a 
temporally varying distance l(t): 

H = p- + ^+V 1 (x 1 -l(t)/2) 
Ira Ira 

+ V 2 (x 2 + l(t)/2) + 2Trhaj tI a s S{x 1 - x 2 ) . (3) 

In this model the ground and first-excited states of 
each trapped atom form a single qubit computational 
basis. We assume distinguishable particles so no 
(anti)symmetrization is required. The trapped particle 
interaction is modelled as an s-wave scattering compo- 
nent of a Van der Walls interaction with scattering length 
a s that reduces upon integration over the transverse de- 
grees of freedom to the above ID interaction, and lou- is 
the frequency characterizing an harmonic approximation 
for the transverse degrees of freedom To avoid prob- 
lems due to degeneracy, the two potential wells have dif- 
ferent individual eigenvalues. Otherwise, when the cou- 
pling between degenerate qubits is switched on, there will 
be fast oscillations between degenerate states regardless 
of how slow the Hamiltonian changes in time. Such a 
double well with different frequencies can be obtained by 
two pairs of counter propagating lasers with wave num- 
bers k, 2k and a bias electric field E. The lasers' relative 
phase and intensities determine details of the double- well 
potential while the constant field can be tuned so that the 
two minima have the same depth. Varying k by changing 
the angle between the incoming beams while increasing 



E and tuning the overall intensity can have the affect 
of changing the distance between minima in time while 
keeping the different trapping frequencies fixed. In Ap- 
pendix [B] we present additional details of an implemen- 
tation of a non-degenerate double well potential in an 
optical lattice. 

When the two traps are at a distance I apart, E mn (l) 
and |m, 77; I) are the instantaneous eigenvalues and eigen- 
vectors respectively: H\m,n;l) = E mn (l(t))\ra,a;l) . 
Here E mn (l(t)) = e\ n + e 2 n + u mn (l(t)), e} n + el is the en- 
ergy of the two atoms in their non-interacting traps and 
u m n{l) is the interaction energy which depends on the 
distance l(t). The asymptotic eigenstates at t = and 
r, are direct-products of the individual trap eigenstates, 
|m) = and \n) = <!>„: \m,n) = \m,n;lo) = \ra) ® \n). 
Such an initial cigenstate evolves into 

M(!;r,0)Kr i ) = ^fi.(OKn') ) (4) 

m'n' 

at time r, where IA(1; t±, to) is the unitary evolution from 
time to to t\ generated by Hamiltonian @ with a gen- 
eral time-dependent distance l(t) and the sum is over 
a complete Hilbert space. A two-qubit gate is a closed 
path l(t), such that 1(0) = 1(t) = lo and the subspace W 
spanned by {|00), |01), |10), 11)} is restored at time t = r, 
with high fidelity A non-operating gate is one for which 
l(t) = lo- The distance lo must therefore be such that 
the interaction energy for states in W can be neglected, 
i.e., u mn (lo) ~ 0, for < n, m < 1. 

IV. OPTIMIZATION OF THE GATE 

We wish to find a path l(t), such that the restriction 
of U(l; t, 0) of Eq. (0} on the computational subspace W, 
U(1;t,0)\w, is equivalent to the required gate of Eq. JQ). 
We do so by finding an objective functional J[l(t)], whose 
minimum is obtained when U(1,t, 0)\w is equivalent to 
the required phase gate. Given a functional, J [I], find- 
ing an optimum l(t) reduces to well established optimal- 
control functional analysis |l4j . An iterative procedure 
is applied; after solving for a given l(t), the functional 
J [I] is evaluated and a gradient search method is used to 
update l(t) as a new, better, trial function. 

In the case of adiabatic evolution, optimal control 
can be used to enhance the fidelity of the gate, but is 
not essential. In adiabatic evolution, as long as no en- 
ergy crossings are involved, an eigenstate evolves to the 
same eigenstate, and Eq. (0} reduces to U(l; r, 0)|m, 77) = 
e~i fa E mn (l(t))dt^ m ^ n y ^ j s t r i v j a Uy satisfied and 
U(i : T 7 0)\w is equivalent to a P(</>) phase gate with = 
J Q T il(l(t))dt, where £1(1) is the acquired controlled phase 
per unit time at distance I given by: 

n(l) = -- (u 0Q (l) + un(0 - u 01 {l) - itio(O) ■ (5) 

Note that £1(1) is not monotonic and can have a maximal 
value £Im at Im (e.g., see Fig. P). 
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In designing a faster gate optimal control becomes 
essential. Restoring the computational subspace after 
rapid switching of the control parameters is not a trivial 
task. After much iterative work, using carefully chosen 
new objectives and numerical methods as detailed be- 
low, we suggest the following fast-approach scheme: (a) 
Change l(t) as fast as possible from Iq to Im, under the 
constraint that (mn; lu\U(l; t(lM.), 0)\m'n'; Iq) oc 6™^ > 
i.e. temporal eigenstates are restored, (b) Let the atoms 
interact for a time w (/)/Qm- (c) Change l(t) as fast as 
possible from Im back to lo, such that the overall evolu- 
tion is diagonal in the computational basis. This is sim- 
ilar to previous suggestions for the phase gate with two 
differences: the approach is not to the smallest possible 
distance but to the optimized distance for acquiring an 
cntanglcmcnt-cffcctivc phase, and the approach and sep- 
aration are not required to take place adiabatically. The 
fundamental bound for the time required to operate the 
gate is reduced to 4>/^Im plus the time required to evolve 
l(t) from Iq to Im and back. As shown below, the time 
required for the approach and descent is reduced in this 
way by an order of magnitude, while high fidelity is main- 
tained. Before presenting the optimal-control results we 
comment on the choice of objective and the parameteri- 
zation of l(t). 

The choice of an objective out of the family of all equiv- 
alent functionals is crucial. One can use the objective 
Jg[1] = — |Trw(W(7, t, 0)G^)| whose minimal value is ob- 
tained when U(1;t,0)\w = G [l^. However, to use this 
we would have to single out a specific gate G £ P (</>), 
whereas a physical system may evolve more naturally to 
another. Instead, we define a new objective Jj>[l] which 
is minimal for any P(<fi) gate: 

J<f,[l] = Icoocoiciocn - exp(i(j))\ , (6) 

where c mn = c™" and c™? n , is defined in Eq. (@J. Similar 
objectives are given by J<p/F and J^j D where 

l 

F=(l/4) l c ™»! 2 ' ( ? ) 

m,n,k,l— 

and 

l 

D = (l/4) l c «™l 2 < ( g ) 

m,n— 

are fidelity measures. F quantifies the unitarity of the 
evolution reduced to the subspace W. This measure es- 
timates how close the evolution is to any two-qubit gate, 
and 1 — F quantifies leakage outside the computational 
basis. The closer the second fidelity measure D is to 
unity, the closer the eigenstates within W evolve to them- 
selves at time r. A stationary path for is necessarily 
a stationary maximal path for F, D, but not vice versa. 

To solve the dynamics with Hamiltonian @, we ex- 
pand a general solution ((t) in the basis of eigenstates 



^mn of the non-perturbed Hamiltonian Hq (a s = 0): 
C(t) = Yl fmn(m mn cxp (~ {el + el) t) , (9) 

mn ^ ' 

® mn (xi,x 2 , 1) = (x! - $1 (x 2 + . (10) 
The Schrodinger equation reduces to 

mf = iif + iMf, (ii) 

where for an operator X, I is short for X = 
cxp (— iHot/h) X exp (iHot/H), and 

U!L=1 J dx^{x-)$l{x + )$l{x-)tf{x + ) ,{12) 

M% n = (* mn \p 2 - Pl \* kl ), (13) 

with pi = —ihd Xi , 7 = 2irfiiOtTa 3 and x± = x ± i. We 
solved the system of differential equations using a stiff 
solver, propagating four vectors of the subspace W. The 
dimension of the matrices was increased until the solution 
converged. The basis we chose is natural for adiabatic 
evolution so larger matrices were required for the non- 
adiabatic simulations. 

The trial functions for l(t) were parameterized by: 

(Q(qi;6i) for < i < ii 
I = I l p for ti « t « t 2 (14) 

l<2(<? 2 ;6> 2 ) for t 2 <t<r 

where 6\ =tjt\ and 8 2 = (r — t)/(r — t 2 ). This describes 
a closed path {l(t) : < t < r}. For t\ « t « i 2 , 
I = l p is fixed; we refer to this as the plateau. The ap- 
proach and departure from the plateau are characterized 
by Q(q]x) = lo + (lo — l p )x 2 (2x—3)+x 2 (x—l) 2 q(x), where 
qi{x) are arbitrary polynomials. With this choice, both I 
and / are continuous. The coefficients of q%, q 2 along with 
tij (and sometimes r) are adjusted by the optimal 

control scheme. This parametrization is suitable for both 
the adiabatic gate and the fast-approach gate. 

V. NUMERICAL RESULTS 

In our numerical example the individual potential wells 
i = 1 , 2 are harmonic with frequencies u>i . The interac- 
tion strength is then characterized by e = 2y / 7ra s Sii;t r / A 
where A 2 = + X 2 and A; = y/ h/muji is the har- 
monic oscillator length for well i. (We only consider 
I > A to avoid nonperturbative effects.) Energy was 
taken in units of fkoi, length in units of A, and time in 
units of l/u)\. The interaction strength was taken to be 
e = 0.05 fvjJi, corresponding to cJtr ~ 10 Wi, a s f=a 0.005 A, 
and A2 = 1.2 Ai. All these numbers were chosen in the 
range of recent experiments in optical lattices. (E.g., for 
87 Rb, a s = 0.005A corresponds to u w 24 Hz). Figure [I] 
shows r2(Z) , the controlled phase acquired per unit time 
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defined in Eq. J5J, where in leading order perturbation 
theory u mn (l) = U™£. £1(1) is non-monotonic with local 
extrema at Mi and M2. Most of the entanglement phase 
is accumulated in our scheme during the plateau. Thus, 
despite the phase lost on the way because the change in 
sign, we expect the best plateau to be at M2 where the 
control-phase acquired per unit time is maximal. 
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FIG. 3: Absolute values of the coefficients f%„(t) (Eq.Q) 
for the optimal adiabatic gate. The upper indexes represent 
the four initial condition /JJ B (0) = Si m Sj n , < i, j < 1. To a 
good approximation, eigenstates are self-evolved at all times, 
as expected in adiabatic evolution. 



FIG. 1: The acquired controlled phase per unit time fi, as 
a function of distance /. The maximal value is at M2 where 
l/X = 1.25; a time of 128/wi is required to obtain a control- 
phase of value 7r. 
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FIG. 4: Absolute values of the coefficients fmn(t)> as in Fig. 3, 
for the fast-approach gate. Eigenstates in W are restored at 
the end of each fast step. 



FIG. 2: The optimal path l(t) obtained after optimization for 
an adiabatic (dashed curve) and fast-approach (solid curve) 
gate (I also plotted for fast-approach). 

We first applied the optimization procedure to an adia- 
batic gate. Good fidelity is expected in an adiabatic evo- 
lution. However, propagating adiabatically, one has an 
error of the order of I. Using the optimization scheme, the 
gate was improved and the error reduced. For the initial 
I trial function we took l p = Im — 1.25 A, and r = 300/wi, 
t\ = t — 1,2 = 100/wi, q — 0. This choice was constrained 



by the requirement that the evolution be adiabatic, i.e., 
that \IM\ <C 1. We expect some phase to be acquired out- 
side the plateau, so we took the plateau time to be smaller 
than 7r/|f2 M | w 128/wi. The maximal value of \l\ is 
0.07A/wi, consistent with adiabaticity. The optimized 
parameters are l° pt = 1.32A = 1.06Z M , *i Pt = 100.47wi, 
t2 Pt = 200.94wi. The optimized I is shown in Fig.0 We 
ended the optimization with J = 0.00003, F = 0.99999, 
and D — 0.99999. Parameter variations of the order of 
10% affect D and F to about one part in 10 5 , while the 
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acquired phase (here <f> = 3.1416) is very sensitive to any 
change in I, in particular 5<f)/4> ~ (St 2 — Sti)/(t2 — ti), 
since cj> ~ ft(l p )(t2 — t\). The coefficients of the evolved 
computational basis are shown in Fig. 3. 

The fast- approach gate is considered next. We ex- 
pect the optimal path for this scheme to be obtained 
when at the end of each fast step (approach or depar- 
ture) all eigenstates in W are restored to themselves 
while the required phase difference is obtained at the 
plateau. To select a trial function we therefore choose 
{h - h) = 128/wi, l M = 1.25A, t x = 11/ u x , r = 150/wi. 
This gives non-adiabatic evolution, max(|Z|) « 0.64 X/u)±. 
The optimized I is the solid curve of Fig. The op- 
timized parameters are l° pt = 1.36 A, t opt = 9.38wi, 

t opt = 140.49 wi, r opt = 150.22 wi. Optimization it- 
erations were stopped when J = 0.006, D = 0.997, 
F = 0.999, and = 3.14. The accuracy for F and D 
is maintained as long as parameters for / are ~ 1% from 
the optimized values. With such variations, the accuracy 
for (j> is 0.01. Varying the plateau time linearly from 
to t opt — t opt while keeping the approach and departure 
shapes fixed, produces a linear curve for <p(t) with slope 
fi(Zp) = 0.0237wi with standard-deviation sw 10~ 4 . With 
this set of plateau times the fidelity measures demon- 
strate oscillatory behavior about their optimized values 
with deviations of about 10 -3 . The coefficients of the 
evolved computational basis arc shown in Fig. 4. Initial 
temporal eigenstates were recovered at the plateau. The 
approach and departure times are reduced by a factor of 
10, leakage outside the computational basis is restricted 
to less than 0.1%, while excellent agreement with an ideal 
P(7r) gate is maintained. 



VI. A FURTHER APPLICATION OF THE 
FAST-APPROACH PROCEDURE 

The fast-approach scheme suggested by our analy- 
sis is essentially model independent. For example, let 
us briefly consider the implementation of such a fast- 
approach scheme for high-fid elity transport in the ex- 
perimental system of Rcfs. pLlj l where the computa- 
tional basis was taken to be the hypcrfine states of 87 Rb 
atoms trapped in an optical lattice in the Mott insulat- 
ing regime. The effective potentials experienced by the 
atoms depended on their internal hypcrfine states and the 
laser polarization. The gate was not implemented adia- 
batically. For approach distances of half a lattice spacing 
and approach times, r, longer than ~ 2-7r/u>, where lo is 
the harmonic frequency at the minimum of the optical 
lattice, the degree of infidelity obtained was ~ 5%. In a 
preliminary study, we applied optimal control to the ap- 
proach step of this gate. The spatial wave packet of the 
atoms (f>(x — R(t), t) can be evolved upon moving the trap 
minima via a path l{t) such that (R— I) + u> 2 (R — I) = —I, 
and the probability to escape the ground state at time t is 
given by P(t) = 1- \{<j>(x- l(t), t)\<f>(x- R(t), t)}\ 2 = 1- 



e ( v (*)/ 4A ) j where R(t) is the center of the wave packet 
andy(t) = R(t)-l(t) with initial condition y(0) = [Tt| . 
Taking \/y 2 (r) + y 2 (r) as the objective, we were able to 
find a smooth path l(t) that reduces the approach time 
by 65%, with y 2 (r) << 10 -12 so that losses were reduced 
from several percent to P(t) ps 2y(r) 2 /4A 2 f» 10 -12 . 
More work is planned to minimize the run time, and im- 
plement the complete gate, including the time interval of 
interaction, yet clearly our scheme is able to significantly 
increase fidelity while somewhat reducing the required 
time for the approach. 



VII. CONCLUSION 



In summary, using optimal-control, we improved the 
adiabatic phase gate, designed a new fast-approach phase 
gate, and demonstrated the new fast-approach scheme for 
a simple model and a recent experiment. More sophis- 
ticated models with additional degrees of freedom could 
further exploit optimal-control techniques. 

A few last remarks are in order. The ID model con- 
sidered here is obtained in reality from a physical 3D 
trap assuming that the ground state of the transverse 
degrees of freedom is maintained through the process. 
A 3D model, where the moving potential remains ID, 
yet transverse excitations are allowed, could be better 
for designing a fast gate. In addition, for any scheme 
to be practical, it is essential that the accuracy of the 
gate depend weakly on small changes in I. Here we saw 
that different fidelity measures are quite robust, while the 
phase is very sensitive to changes in I. As a result the 
gates can be easily adjusted to give a P(4>) gate with any 
4> by changing the time duration of the plateau. In prac- 
tice, feedback learning techniques may be used to assure 
stability and to adjust the phase. Finally, for a complete 
scheme of computation, one also needs single qubit gates. 
This is another point at which degeneracy causes prob- 
lems, as the computational basis of the two oscillators 
needs to be un-entangled. While these points deserve 
further study, they were disregarded here, as our main 
focus was that one could use optimal-control techniques 
and a fast-approach scheme for the implementation of 
the control-phase gate in models of cold atoms in an op- 
tical lattice, significantly reducing the gates time while 
maintaining high fidelity. 
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APPENDIX A: ENTANGLEMENT POWER OF 
THE CP(</>) GATE 

An important property of the CP((j>) gate (and hence of 
any P(4>) gate) is that the phase <fi can be used to measure 
the entanglement power of the gate. To see this consider 
a general two qubit state \^f) = J2 m n=o a mn\fn, n). The 
action of CP ((f) on this state is, 



A=( a 00 a 01 

~~ V a io an 



A = 



a 0Q 0-01 

a 10 cran 



(Al) 



where a = cxp(ic/>). Generally if p is the reduced density 
matrix of a bipartite pure system, the usual entangle- 
ment measures of p are given by the purity (linear en- 
tanglement) measure Mi = — Trp 2 , or the Von Neumann 
entropy M 2 = — Tip log 2 p. There are important mathe- 
matical features of interest when one chooses a measure, 
such as convexity. Here we are only interested in a mea- 
sure of the degree of entanglement. The reduced density 
matrix of two qubits can be written as 



p = P T + (l-p)(I-V), 



(A2) 



where < p < 1, V is the one dimensional projection 
onto the p eigenstatc of p, and / is the identity in the 
two level system sub-space. For such a state 

-[p 2 + (1 - pf] , (A3) 
-plog 2 p-(l-p)log 2 (l-p). (A4) 



Mi 
M 2 



The line ordering imposed by these measures is the same 
line ordering obtained from the determinant of the re- 
duced density matrix: det(p) = p(p — 1). We therefore 
use the change in this determinant as a measure for the 
entanglement induced by the gate: 



A(*;0) = det(Ii+) - det(AA t ) (A5) 
= 2!ft (a Q0 an a 10 ooi (1 - a)) (A6) 

= m (-i rexp (i( 7 + ^ sin(^)J (A7) 

where rexp(^7) = aoo an «io «oi- The maximal change 
is obtained when r and sin(7 + </f>/2) are maximal, i.e., 
when |a mn | = 1/2 and 7 = (tt — (f>)/2. The entangling 
power of the CP(</>) gate and hence of any P(</>) gate is 



max«(A(tt,0)) = jsin(|; 



(A8) 



Maximum entanglement power is obtained for <fi = tt 
(also for the CN gate). 



APPENDIX B: ASYMMETRIC DOUBLE- WELL 

Two counter-propagating laser beams at a fundamen- 
tal frequency and at its second harmonic, with wave num- 
bers k and 2k, and intensities, produce an effective stand- 
ing wave optical potential of the form 

V{k,x) = V [cos 2 (kx + 5) + acos 2 {2kx)], (Bl) 



where Vq is proportional to the intensity of the funda- 
mental, a is determined by beam intensities and detun- 
ing ratios, and S is determined by the relative phase be- 
tween the fundamental and second harmonic. The wave 
number k can be tuned by manipulating the angle be- 
tween the right and left laser beams. The effective po- 
tential thus formed is an array of double wells, with a 
determining the relative heights of the well minima. For 
5 > 0, the wells are non-symmetric. Let Xi(k), x 2 (fc) 
denote the positions of the minima of a specific double 
well, where x\(k) < a; 2 (/c). The greater 5, the larger 
V(x2{k)) relative to V(xi(k)), as well as V"(x2{k)) rela- 
tive to V" (xi(k)). One can add a linear potential of the 
form 



V{x,k) = -A(k)x 



(B2) 



to V(k, x) so the resulting potential, V e s , will have equal 
well minima, V b q (a; 2 (/c)) = V c g(xi(k)). This requires 



A(k) := Vc 



V(x 2 (k))-V(x 1 (k)) 
x 2 (k) ~ xi(k) 



(B3) 



As a result of the additional linear potential, there is an 
increase in the ratio of the second derivatives at the well 
minima of V e g. Note that the original potential V(k,x) 
is a function of kx, V(k,x) = f(kx), thus 



V(k',x) 

and hence 

A(k') 



f(kjx) = V(k,jx), 
kxi (Jz) , 



V(k',x 2 (k'))-V(k',x 1 (k')) 

k' (x 2 (k') - Xl (k')) 
V(k,x 2 (k))-V(k, Xl (k)) 



A(k) 



k (x 2 (k) - xi(k)) 



(B4) 
(B5) 

(B6) 
(B7) 



A(k)/k is independent of k and l|B2|l can be written as 
V{k',x) = -V Ck'x (B8) 



where C is a constant that may be determined by l|B3|) 
for any k. Thus, 

V eS (k,x) = V [cos 2 (kx + S) + acos 2 (2kx) - Ckx] ,(B9) 

is a function of kx. For such a potential, variations of 
k scale x. The distance l(k) between the two minima in 
a double well is thus a parameter scaled by k, l(k') = 
^-Z(fc), while keeping the values of potential values at 
minimal points unaltered as well as the ratios of second 
derivatives at the minima. To maintain the values of 
second derivatives at the minimum constant as one varies 
k, one can modify the laser intensities to depend on k 
such that 



V (k) :=|1. 



(BIO) 
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